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Abstract. We calculate a nondissipative spin current and show that it can flow with or 
without a charge current. We consider a two-band model which can be applied to the 
description of Fe-based pnictides in coexistence regime of superconductivity and spin- 
density wave. Using quasiclassical Green’s functions approach and tunneling Hamiltonian 
method we show that there exists a possibility to switch off the Josephson current while 
leaving the spin current finite. Moreover, it is possible to have the critical Josephson current 
and the critical spin current being proportional to each other, thus giving a possibility to 
measure the spin current via the Josephson current. The underlying mechanism is the 
interfering hopping of electron-hole pairs between different bands of the superconductors 
composing the junction. This is an intrinsic property of the system and provides a unique 
and natural way to utilize junctions made solely of pnictides in promising applications in 
spintronics devices. 
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1. Introduction 

Control of spin currents without a net charge transfer 
constitutes a novel possibility to transfer and process in¬ 
formation, and study questions of fundamental physics, 
and also offers a promising variety of applications [1, 2]. 
Advantages, aside from faster switching times and lower 
power consumption as compared to conventional de¬ 
vices [3], are a higher integration density [4] and a high- 
bandwidth information transfer [2]. 

A major challenge, though, are creation, detection, 
and control of spin currents. There are various possi¬ 
bilities discussed in the scientific community including 
the use of circularly polarized light and electrical spin in¬ 
jection from magnetic contacts [1, 2], or creation of spin 
current in ferromagnet-semiconductor junctions [5] and 
its enhancement in ferromagnetic insulators [6]. Since, 
for practical implementations, generation and creation 
of spin current should not require strong magnetic fields 
and interfaces between semiconductors and ferromag- 
nets, Sharma proposed a method of pure spin pumping to 
create spin current using a two-dimensional electron gas 
and a quantum cavity [3]. However, the detection of spin 
currents with arbitrary polarization direction, still consti¬ 
tutes a nontrivial task. 

Suggestions to produce, measure, and control pure 
spin current in devices mostly based on semiconduc¬ 
tors include theoretical studies of the influence of elec¬ 
tric fields [7, 8], the optical quantum interference injec¬ 
tion [9], and mechanical torque observation [10]. More¬ 
over, it has been proposed that spin current may arise in 
mesoscopic hybrid structures due to interfacial spin-orbit 
scattering [11] or at twin boundaries of noncentrosym- 
metric superconductors [12]. Spin-orbit coupling is an 
important ingredient also in a proposal to use a single¬ 
channel voltage probe to measure spin current in lateral 
heterostructures [13]. 

Aside from coupled ferromagnets [14], spin current 
and the Josephson effect have been investigated in triplet 
superconductor junctions with an insulating [15] orferro- 
magnetic [16] layer. An interesting interplay between fer¬ 
romagnetism and superconductivity influencing charge 
and spin current has been found in a structure consisting 
of two nonunitary ferromagnetic spin-triplet supercon¬ 
ductors separated by a thin insulating layer [17]. More re¬ 
cently, spin current has been investigated in heterostruc¬ 
tures consisting of triplet and singlet superconductors in 
contact with a ferromagnet [18], or in a Josephson junc¬ 
tion with double layer ferromagnets in contact with s- 


wave superconductors [19] where, by analogy to the long- 
range component of charge Josephson current [20], a 
long-range spin current may arise as a consequence of the 
presence of spin-triplet Cooper pairs formed by electrons 
of equal spin induced by the long-range proximity effect 
inside the ferromagnet [21]. A similar phenomenon, i.e., 
the long-range penetration of the antiferromagnetic or¬ 
der parameter (and, correspondingly, of the spin current) 
into a ferromagnet, may occur in a purely magnetic struc¬ 
ture [22]. 

Exceptional sources of spin currents might be pro¬ 
vided by single-molecule magnets [23] and topological 
superconductors being nowadays under intense investi¬ 
gation [24]. 

However, the discovery of superconductivity in 
so-called iron based pnictides [25] revealed a novel 
class of high- T c superconductors with a multiple band 
structure where superconductivity may coexist with 
spin-density wave, thus constituting a system which is 
intrinsically superconducting and magnetic at the same 
instance. This feature has great potential for fundamental 
research as well as for devices based on involvement of 
superconducting and spin degrees of freedom. 

In this Paper, we study the spin current in a Joseph¬ 
son junction composed of two-band superconductors 
with a spin-density wave separated by a thin insulating 
layer. We use a model which is widely accepted to de¬ 
scribe qualitatively the coexistence region of supercon¬ 
ductivity and spin-density wave in Fe-pnictides forming 
a Josephson junction [26, 27, 28, 29, 30]. Novel to inves¬ 
tigations on spin current is the interplay of superconduc¬ 
tivity and magnetism and also the possibility of tunnel¬ 
ing of electrons and holes (and pairs of them) between 
nonequal bands. Despite the relative simplicity of the 
model, we find that, aside from the usual Josephson cur¬ 
rent, a nondissipative spin current can flow through the 
junction. Moreover, choosing proper conditions one can 
separate charge and spin degrees of freedom both intrin¬ 
sic to the considered system. After the introduction of the 
model we provide expressions for the spin current and the 
corresponding expressions for the Josephson current in 
the case of symmetrical and asymmetrical junctions dis¬ 
tinguished by the pairing properties of materials on each 
side of the junction (see next section for details). Finally, 
the obtained results are analyzed and an experimental 
setup is proposed to test them. 
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2. Model 

We consider a two-band model for pnictides [26, 27, 
28, 29, 30] suitable to qualitatively describe coexistence 
region of superconductivity and spin-density wave in 
these materials. Introducing appropriate operators [30], 
the Hamiltonian can be written as 

( 1 ) 

2 p 

with the 8 x 8 matrix H = H^ n + Hf c + H s dw consisting 
of the corresponding parts, i.e., l\ lt) - c, r (p)Xo .30 - pA/son, 
Hg C - £(p)A'X(o, 3 )i 3 - A"X(3,o) 23, ffsdw = tnX n 3 , where the 
order parameter A = A' + iA" is related to the supercon¬ 
ducting energy gap Ai in the hole pocket as A = A* 
and the cases of s++- and s+_-pairing are distinguished, 
i.e., in the case of s++-pairing, we have for the SC or¬ 
der parameter in the electron band A 2 = Ai, while for 
s+_-pairing, A 2 = — Ai (the magnitudes are taken to be 
equal). Moreover, we introduced the matrices X mna 
denoting the Kronecker product of the corresponding 
Pauli matrices in the band, particle-hole, and spin 
space, respectively, X mna - p m T n <j a [30, 31]. The disper¬ 
sion relations can be linearized near the Fermi energy, 
1 , 2 (p) = +c(p) + p with (f(p) = vi p, the Fermi velocity v F , 
and the nesting parameter p-po + pc/, cos ( 2 (p), where po 
describes the relative size difference of electron and hole 
pockets and p$ controls the ellipticity of the electron 
pocket. Moreover, the spin-density wave order parame¬ 
ter is defined as m — ra q= 0 + m* =0 , see [26, 27, 28, 30] for 
further details. 

In terms of the introduced operators one de¬ 
fines the Green’s functions, e.g., the retarded and 
Keldysh Green’s functions, and, correspondingly, the 
quasiclassical Green’s functions for each component, 
g = ^ /d<f(Xo 3 oG). The thus introduced quasiclassical 
Green’s functions obey the generalized Eilenberger equa¬ 
tion [30], 

v F Vg+[w„Xo 30 + iA,g] =0, (2) 

supplemented with the normalization condition, g 2 = 1 . 
Here, the matrix A = (- pXxw + H^ c + H s dw) A 030 , and 
(i) n - (2 n + 1) n T are the Matsubara frequencies. 

The solution g ,0) = g R in the case of real supercon¬ 
ducting gap and the magnetization of the spin-density 
wave being oriented along the z axis has the form 

g|°2 = g030 + glOO + gl23 + g213 + g300 + g323 for the S+-- 

pairing, and g™ = | 023 +1030 +1123 +1130 +1213 +1300 
for the s ++ -pairing. If in contact, the phase difference 
of the SC condensates, cp, and the mutual orientation of 
the magnetization of the SDW in the leads, a, can be 
accounted for by the functions gi( r j, in the left (right) 
lead, respectively, which are related to g t0) via the 
unitary transformation gi (r) = R la S ±ll) g <(]> h\ a S j ±(fl , where 
S+cp - exp(±iX 33 o<p/4) and R ±a = exp(+iX 33 ia/4). They 
can be called the rotation matrices in Gor’kov-Nambu 
and spin spaces, respectively. 
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Figure 1. (Color online) Considered setup. The phases in the 
superconductors are assumed to be ±cp/2, and the angles between the 
z axis and the magnetization vectors of the SDW are ±al2. 



Figure 2. (Color online) Schematic tunneling through an interface 
between two two-band materials. 

To calculate the spin current we employ the tunnel¬ 
ing Hamiltonian method, i.e., we consider a tunnel junc¬ 
tion (see figure 1 ) composed of two two-band supercon¬ 
ductors with a spin-density wave where each SC is de¬ 
scribed by a Hamiltonian of the form (1), while transitions 
of electrons between them are accounted for by the tun¬ 
neling Hamiltonian [22, 32], 

^, = i£[C r t .rCi + h.c.], (3) 

2 p 

where the transfer matrix can be written as 3~ — U~q C Xo 30 
with^Fjjc = -5 + X3oo-(5AXooo+i[[R(<5l2)X210 + 3("5l2)A22o]> 
where 3~± — [ST\\ ±322)12 with tunneling amplitudes STu 
between equal bands, and ,Tj 2 describes tunneling be¬ 
tween bands 1 and 2, cf. figure 2. Without restriction, 
the tunneling elements between equal bands can be sup¬ 
posed as real quantities (we assume that the matrix ele¬ 
ments 3~ij do not depend on momentum p), but the ma¬ 
trix element 3\ 2 is, in general, complex [32]. 

The tunneling current in equilibrium is the nondis- 
sipative Josephson current 7j. It is given by 

h = i (2n T) £ Tr ( X 330 [ .T qc g\ lt»).t qc , g r M ]), (4) 

Ik (J) 

where vi r and g] , (w) are, correspondingly, the densities 
of states and the quasiclassical Green’s functions in the 
Matsubara representation in the left and right leads, 
respectively, the angle brackets mean averaging over all 
directions of p, ‘Tr’ denotes the operation of taking trace 
of a matrix, and e is the elementary charge. 
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In analogy to the Josephson current, the spin current 
can be written as 

h P = i n ^ V ' Vr (2n 7) £Tr(X 33 i [.%cg\ M-Tqc, g r (co )]), (5) 
IV (1) 

where p b is the Bohr’s magneton. 

3. Spin current 

We calculate the spin current (5) for three possible ar¬ 
rangements, i.e., the symmetric s ++ Is ++ , respectively, 
s + -Is+- junctions, and an asymmetric s++ 7s + _ contact 
(here, 7’ denotes an insulating layer and ‘s’ —a supercon¬ 
ductor with corresponding pairing symmetry in coexis¬ 
tence regime with a spin-density wave). The Josephson 
current (4) has been calculated elsewhere [32], but for 
convenience the corresponding expressions are provided 
here. 

(i) Spin current. In both cases we obtain the Josephson- 
like expression, 

7 S p = 7 C sin a , (6) 

with the critical current having a different form 
depending on the configuration. 

For the symmetric setup it reads 

Ic = lo + Icp cos (f> (7) 

with 

7 0 oc mim r [51i522±5R(^jl)]MBV 2 , (8) 

70 cx mi/7t r AiA r [^ 11^22 ± \$ 12 \ 2 ] pbv 2 > (9) 

where the upper sign denotes the s ++ Is ++ junction 
and the lower—the s+_7s+_ junction. 

For the asymmetric s++7s+_ setup, the critical spin 
current is proportional to 

7 c oc mim r [51i522 + J7(^j 2 2 )]pBV 2 (10) 

and does not depend on the phase difference 
between the superconductors. 

Similar expressions for the spin current were ob¬ 
tained in [33] for a FFLO-like superconductor with 
a helimagnetic molecular held. In contrast to the as¬ 
sumptions of the authors of [33], in our case, a coex¬ 
istence of superconductivity and magnetic order is 
intrinsic to the considered model and confirmed by 
experiments on Fe-based pnictides. 

(ii) Josephson current. In the case of a symmetric 
s ++ Is ++ , respectively, s+_ Is+- junction, we ob- 


tained [32] the Josephson relation 

7j = 7 CiJ sin0, 

where 

(11) 

7cj = 7 0lJ + 7 a j cos a 

(12) 

depends on the mutual orientation of the magneti¬ 
zation in the spin-density wave in the leads, with 

7 0 ,j cx Ai A r [3\\ + ,3T 2 2 2 + 23R(,7j|) ] ev 2 , 

(13) 

7a,j cx mim r AiA r [5'n^2±l^i2l 2 ]ev 2 . 

(14) 


Again, the upper sign represents the result for the 
s ++ Is ++ junction and the lower—for the s+_7s+_ 
junction. 

In the case of an asymmetric setup, an unusual so- 
called 0-junction can be realized [32], 

7 c ,j = I\ sin0 + / 2 cos0, (15) 

with Ii oc AiA r (77j 2 an< ^ ^ AiA r 3(,^|). 

4. Analysis 

First of all, as can be seen from the definition of 
corresponding coefficients in the expressions for spin 
currents (8) and (9), switching off superconductivity (in 
pnictides, this can be done if considering the sample 
at a temperature above the superconducting transition 
temperature but below the Neel point), we obtain a result 
resembling the situation considered in [22] and [14], 
i.e., a nondissipative spin current flowing through the 
junction due to a misalignment of the magnetization 
directions in the banks, where the spin current is carried 
by quasiparticles of same type (electrons or holes in the 
first case and electrons in the second case). In our model, 
there appears an additional component (proportiopnal 
to ST\ 2 ) of spin (and also charge) current due to presence 
of two bands and due to possibility of hopping between 
different bands, i.e., electron-hole pairs hopping occurs 
with a total spin of 1 and no total charge transfer. 

However, in coexistence regime, in the symmetric 
junctions, there appear corresponding contributions to 
the spin and charge currents due to pair hopping of 
Cooper pairs and Cooper-like electron-hole pairs carry¬ 
ing spin but no charge, leading to nontrivial mixture of 
charge and spin degrees of freedom. These contributions 
are represented by the phase dependent part (9) propor¬ 
tional to A[A r for the spin current, and by the angle de¬ 
pendent part (14) proportional to m\m r for the Josephson 
charge current. 

In this situation, the most striking result is the pos¬ 
sibility to separate the charge and spin degrees of free¬ 
dom in all considered combinations. Notably, this is 
an intrinsic property of the considered physical system, 
as opposed, e.g., to [19], where the authors investigated 
a Josephson junction made of two conventional super¬ 
conductors with two diffusive ferromagnetic layers sand¬ 
wiched between those, by means of quasiclassical Green’s 
functions and solving the Usadel equation for these im¬ 
posing appropriate boundary conditions. In the present 
case, coexistence of superconductivity and magnetism is 
intrinsic to the considered system, but in spite of that, 
they can be separated, as we now show based on an anal¬ 
ysis of corresponding coefficients of spin (7o, /,/, and 7 C ) 
and charge (7oj, la,] and 7 c j) currents, see figure 3. 

In the case of an asymmetric s ++ Is+- junction, 
the Josephson current vanishes if 3~ xx - TT, 2 = 0 and 
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Figure 3. (Color online) Spin [solid red line denotes I c and Iq 
(respectively I ^)] and Josephson [short-dashed black (/ij) or long- 
dashed blue (Iqj) and dash-dotted green [I a j) line, respectively] 
currents on the normalized tunneling coefficient r = 3 22 /^"ll in the 
asymmetric (a) and symmetric (b) case in arbitrary units. Other 
parameters are $R(#i2/^il) = 0 and = 1. For r = — 1 the 

Josephson current vanishes, whereas the spin current stays finite. Note 
also the possibility of a sign change for the spin current as well as for the 
Josephson current. 


3(5^1) = 0. The spin current stays finite provided 

Sii&n *-H(#ii). 

In the case of a symmetric s++7s++, respectively, 
s+-Is+- junction, the Josephson current vanishes if 
3~\ 1 ± \3~n 1 2 = 0 and 3~^_ + 3~£ 2 ±25R(/Jj|) = 0. In con¬ 
trast, the spin current stays finite provided f 0, 

since even if the I,/, part of the spin current vanishes due 
to the first condition here, the independent of <p part stays 
finite. 

Moreover, in the symmetric case it is possible 
to make the Josephson current vanish even without 
controlling the tunneling elements. Choosing the angle a 
of the mutual orientation of the magnetization directions 
in the spin-density wave to be a = arccos( + 7oj/7 a j), 
provided 7oj < 7 a j, which is possible, e.g., if m is 
sufficiently large, the Josephson current vanishes, while 
the spin current stays finite. 

Nevertheless, there exists another effect based on 
controlling the tunneling elements in the symmetric 
junctions. As is evident from the expressions of Iq (8) 
and 7o,j (13), if the relation 3\\ = ±^22 holds for s++Is ++ , 
respectively, s+-Is+- junction, the spin and Josephson 
currents are related with each other, 

7 sp = 7<^ cos (/> sin a, (16) 

7j = 7 a j cosasincji, (17) 

via to the relation /,/, = ®-7 a j which follows from the 
comparison of the corresponding coefficients. This 
offers a possibility to measure the spin current via the 
Josephson critical current, i.e., measuring the value of 7j, 
the phase (p , and the angle a, one obtains 7 sp from this 
relation. 

5. Experiment proposal 

As the latter result is the most plausible from the 
experimental point of view, it may be tested by means 
of a nonlocal spin current measurement using a device 
adopted from [34, 35] based on the spin Hall effect. Here, 
the spin current is injected replacing the ferromagnetic 
electrode used to inject spins in figure 1 of [34] by 



Figure 4. (Color online) Suggested experimental setting to detect the 
spin current in the considered Josephson junction. The Hall cross (HC) 
consists of a two-band material; the red electrode (FM) is a ferromagnet 
used for sample characterization [34]; the spin current I sp is injected 
via the contact of the considered S s ^ w / S s d w junction (the two banks are 
shown in green and blue) with the Hall cross. 

the considered Josephson junction, see figure 4. The 
used materials should be accordingly adapted, e.g., the 
aluminium Hall cross should be substituted with a Hall 
cross made of corresponding two-band material with a 
hole and an electron band to capture the spin current 
promoted by the electron-hole pairs with zero total 
charge. 

6. Conclusion 

We considered a two-band model for Fe-based pnic- 
tides and investigated, based on quasiclassical Green’s 
functions approach and tunneling Hamiltonian method, 
Josephson junctions of two types, i.e., a symmetric con¬ 
figuration of superconductors with the same pairing sym¬ 
metry (s ++ 7s++ and s + _7s+_), as well as with different 
pairing symmetry (s++7s+_). We found that it is possi¬ 
ble to separate the charge and spin degrees of freedom 
by having a finite spin current and vanishing charge cur¬ 
rent at the same time. The reverse effect is also possi¬ 
ble but it is not discussed in details. Important ingredi¬ 
ents are the coexistence of superconductivity and spin- 
density wave and the hopping between nonequal bands. 
Note that the possibility to obtain a superconductor with 
s ++ pairing symmetry out of s+_ pairing by adding impu¬ 
rities has been investigated in [36]. The discovered phe¬ 
nomenon may provide a reliable tool for control of spin 
currents, especially, as a detector or a source of spin cur¬ 
rents based on measurement of the Josephson current, 
thus enabling its application in calibration of correspond¬ 
ing devices. Finally, we suggested an experimental setup 
which is suitable test the obtained results via a nonlocal 
spin current measurement. 
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